^APPLICATION NO. 10/806,016 
^^'nv^^y INVENTION: Multi-scale code division frequency /wavelet multiple 



access 

INVENTORS: Urbain Alfred von der Embse 



Clean version of how the SPECIFICATION 
will read. 



1 



APPLICATION NO. 10/806,016 

INVENTION: Multi-scale code division frequency/ wavelet multiple 
access 

INVENTORS: Urbain Alfred von der Eitibse 

5 

This patent apolicaiton is a continuation in part of 
application 09/826,118 filed on 01/09/2001, application 
09/526.117 filed on 01/09/2001, and application 10/266,257 filed 
on 10/08/2002. 

10 

BACKGROUND OF THE INVENTION 

I. Field of Invention 

The present invention relates to both orthogonal frequency 

15 ■ division multiple access OFDMA, to orthogonal Wavelet division 

multiple access OWDMA, to code division multiple access CDMA, and 
to multi-scale code division multiple access MS-CDMA, for 
cellular telephone and wireless data communications with data 
rates up to multiple Tl (1.544 Mbps), El (2.048 Mbps) , Sonet, 

20 Ethernet, and higher (>10 Gbps), and to optical CDMA and optical 
OWDMA. Applications are to wire, wireless local area, wide area, 
mobile, point-to-point, and satellite communication networks. 
More specifically the present invention relates to a new and 
novel means for combining MS-CDMA with OFDMA, to a new and novel 

25 OWDMA which is an orthogonal multi-resolution complex Wavelet 

multiple access generalization of OFDMA, and to a new and novel 
means for combining MS-CDMA with OWDMA. This new architecture 
MS-CDMA OFDMA/OWDMA is an attractive candidate to replace current 
and future OFDMA applications and CDMA applications . 

30 

II. Description of Related Art 

Current OFDMA art is represented by the applications to the 
wireless cellular communications standards IEEE 802.11a, IEEE 
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802. llg, IEEE 802.15.3a, IEEE 802.16. OFDMA uses the Fourier 
transform basis vectors as the orthogonal channelization vectors 
for communications with each basis vector multiplied by a symbol 
which is encoded with a data or pilot signal word. 

5 

The discrete Fourier transform DFT implemented as the fast 
Fourier transform FFT is defined in equations (1) . Step 1 defines 
the digital sartpling interval T over time, the sampling instants 
t=iT where i is the time index and where the sampling time 1/T is 
10 at least equal to the complex Nyquist sampling rate to prevent 
spectral foldover. Step 2 is the FFT of the complex baseband 
transmitted signal z(i) for the data block 

and step 3 defines the NxN orthogonal complex DFT matrix E row 
vectors E(k) which are the DFT harmonic vectors or basis vectors 
15 or code vectors or channelization vectors. Step 4 defines z{i) 
for one data block and is equal to the inverse FFT transform 
FFT'^ of the user symbols x(k). 



Unweighted DFT encoding (1) 
20 1 Sanpling interval of DFT 

t -> 0 T .... (N-2)T (N-l)T 

i ^ 0 1 N-1 N-1 

where 1/T ^ complex Nyquist sample rate 
25 2 FFT of z(i) 

X(k) = FFT[ z(i)] 

= Zi E(k,i) z(i) 

= Zi exp(-j27C ki/N) z(i) 



30 3 DFT orthogonal harmonic code matrix E 

E = NxN DFT orthogonal DFT matrix 

= [E(k,i)] matrix of elements E(k,i) 
E(k,i)= exp( -j27C ki/N ) harmonic k, time index i 
E(k) = harmonic k basis (code) vector 
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= [E(k,0), E(k,l), . . . ,E(k,N-l)] 
EE* = N I 

E(k)E*(k' ) = 5(k-k' ) N I 
where I = NxN identify matrix 
5 E*= conplex conjugate transpose of E 

8(k-k')= Dirac delta function 
= 1 for k=k' 
= 0 otherwise 

10 4 Transmitted DFT complex baseband signal z{i) for 

one data block 

d(k) = data modulation for user k 

= encoded amplitude A(k) and phase cp(k) 

x(k) = transmitted symbol encoded with d(k) 
15 = A(k) exp(jcp(k) ) 

z(i) = FFT-^ x(k) ] 

= N"^ Zjc x(k) E*(k,i) 



20 OFDMA for IEEE 802. llg in reference [1] is illustrated in 

FIG. 1- The channelization filter h(f) 1 covers a 20 MHz 
frequency band 2 assigned to OFDMA. Plotted is the power spectral 
density PSD=|h(f)|^ of this channelization filter h(f). A N=64 
point fast Fourier transform FFT covers this band 2. Consistant 

25 with the IEEE specification, FIG. 1 refers to the DFT which is 
identical to the analog Fourier transform FT since it is the 
sampled data version of the FT. It is convenient to consider the 
DFT in this invention disclosure as the digital format for the 
FFT. The DFT frequency spectrum 3 consists of N=64 equally spaced 

30 filters 4 across this 20 MHz band. Filter spacing is equal to the 
DFT output rate 1/NT == 0.3125 MHz = 20MHz/64. The DFT time pulse 
p(t) 5 is NT=3.2|is in length and the total DFT period 6 is 4.0^3 
which allows a O-S^is guard time for p(t). 
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Throughout this invention disclosure it will be understood 
that the FFT fast algorithm will always be used to implement the 
DFT and the inverse FFT"^ fast algorithm will always be used to 
implement the inverse DFT'^. 

OFDMA transmitter encoding of the OFDMA waveform in FIG.l 
is defined in equations (2) . Step 1 lists the parameters and 
definitions and step 2 defines the time domain weighting. Step 3 
is the complex baseband transmitted signal z(i). 
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OFDMA encoding for transmitter (2) 

1 Parameters and definitions from 1 in equation (1) 
and 

h(i) =20 MHz band filter impulse response 
15 p(i) = impulse response of the DFT waveform 

= real weighting function in 6 in FIG. 1 
N =64 point DFT 

1/T =20 MHz sample rate for DFT 
> complex Nyquist rate 
20 NT = 3.2|is DFT length 

1/NT = 0.3125 MHz DFT output rate 
= DFT channel separation 
= DFT tone spacings 
52 channels are used: 4 pilot, 48 data 
25 12 guard band channels for rolloff of the h(k) 

2 Pulse p and band filter h weighting for DFT basis vectors 

£ = p©h, convolution of p and h 

= filter transfer function in time domain 
30 for the combined h and p filters 

= [ e(0), Ed), . . p(N-l) ] 

3 Transmitted OFDMA encoded baseband signal z(i) 
for one data block 
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d(k) 



= data modulation for user k 



x(k) 



= encoded amplitude A(k) and phase (p(k) 
= transmitted symbol encoded with d(k) 
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z(i) 



= A(k) exp(j(p(k) ) 

FFT'^[ e(i) x(k) ] 

N"^ ZKe(i) x(k) E*(k,i) 



OFDMA for IEEE 802. llg has the strict orthogonality of the 
DFT(FFT) replaced by cross-correlations between the 48 channel 
10 tones and other impacts due to the band channelization and pulse 

weighting p©h plus the time errors At and frequency errors Af 
from synchronization errors, multi-path, propagation, and 
terminal stresses. These inpacts on orthogonality are low enough 
to allow OFDMA to support higher values for the symbol signal-to- 

15 noise ratio S/N in the detection band that are required for 
higher order symbol modulations. The highest order symbol 
modulation currently is 64 state quadrature ait^litude modulation 
64-QAM corresponding to 6 bits per symbol where 6=log2(64) and 
log2(o) is the logarithm to the base 2- With rate 3/4 

20 convolutional coding the highest information rate is 4.5 
bits/symbol = 6x3/4. Required S/N at a BER=1.0e-6 is 
approximately S/N- 19 dB. 



25 frequency band at a symbol rate equal to 0.25 MH2=l/4.0 fis from 6 
in FIG. 1 and for a maximum information rate equal to 4.5 
bits/symbol this equals a burst rate of 54 MBps = 4.5x48x0.25. 
Some spread spectrum properties are realized by hopping the 20 
MHz band, shuffling the channel assignments over the 48 available 

30 channels for a user assigned to several channels in order to 

spread his transmissions over the band, and for "flash" ODFMA by 
a random hopping of each user channel across the 48 available 
channels within the band. 



OFDMA for IEEE 802. llg provides 48 channels over a 20 MHz 
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OFDMA receiver decoding of the OFDMA waveform in FIG.l is 
defined in equations (3) for the receiver and derives estimates 
of the transmitted symbols by implementing matched filter 
detection in the receiver. 

5 

OFDMA decoding for receiver (3) 
OFDMA decoding for one data block derives estimates 
x(k) of x(k) from the receiver estimates z(i) of 2(i) 
x(k) = FFT[ z(i)©p ] 

10 = Si z(i)©p_E(k,i) 

= x(k) 

Current CDMA spread spectrum art is illustrated by the 
waveform in FIG. 2 which describes the waveform for full band 
15 CDMA communications over the band B 9 which is the output range 8 
of the band filter h(f) 7. The CDMA chip rate 1/Tc 12 is less 
than the available frequency band B to allow the chip frequency 
spectrum p(f) 10,11 to roll off. As defined 9 in FIG. 2 the band 

B and chip rate are related by equation B=(l+ a) /Tc where a is 
20 the bandwidth expansion factor close to a=0.25 for high 

performance communications. Frequency spectrum p(f) 10 for the 
CDMA communications is essentially equal to therepresentative 
time pulse p(t) is a square-root raised cosine pulse which can be 
used for high performance communications to obtain a reasonably 
25 flat spectrum with a sharp rolloff at the edges to enable the 

chip rate 1/Tc to be reasonably close to the available frequency 
band B. 

Chip rate 1/Tc is the CDMA total symbol rate. The users 
30 could be at different data rates but this and other architectural 
variations do not limit the scope of this invention. Power is 
uniformly spread over the CDMA pulse waveform spectrum p(f ) . 
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It is self evident to anyone skilled in the CDMA 
communications art that these communications mode assumptions 
are both reasonable and representative of the current CDMA art 
and do not limit the applicability of this invention. 

5 

CDMA encoding of the waveform in FIG. 2 for the transmitter 
is defined in equations (4). Steps 1,2 define the CDMA 
transmission and parameters. Step 3 defines the user symbol 
x(u). Step 4 is the set of Walsh orthogonal channelization codes 
10 w(u) and step 5 is the pseudo-random PN covering or spreading 
code. Step 6 defines the complex baseband signal z(t) as the 
waveform £(t-nTc) multiplied by the encoded sum over u and n. 



15 CDMA encoding for transmitter (4) 

1 CDMA transmission 

chip 0 chip 1 chip Nc-1 

t t t 

t 0 Tc (Nc-l)Tc 

20 n -> 0 1 Nc-1 

i 1 1 

2 ^ £(t) £(t-Tc) E(t-(Nc-l)Tc) 



where p(t-nTc) = p(t-nTc) ©h(t) 
25 = convolution of p(t-nTc) and h 

= filter transfer function in time domain 
for the combined p(t-nTc) and h filters 



2 Parameters and definitions 
30 Nc = Number of users and orthogonal code chips 

Tc = CDMA chip length or repetition interval 
1/NcTc = User symbol rate 



3 User complex signal x(i) 
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d(u) = data modulation for user u 

= encoded amplitude A(u) and phase (p(u) 

x(u) = transmitted symbol encoded with d(u) 
= A(u) exp(j(p(u) ) 

4 Walsh orthogonal channelizaation code matrix W 

W = Code matrix, Nc rows of Nc code vectors 

= [ W(k,n) ] matrix of elements C(k,n) 

W(u,n)= +/-1, chip n of vector u 

W(u) = code vector u, row k of W 

= [W(u,0), W{u,l), . . . ,W(u,Nc-l)] 

5 PN covering (spreading) code P(n) for chip n 

P(n) = exp(j(p(n)) 

6 Transmitted CDMA complex baseband signal z(t) 

z(t) = Nc"'Iu In P(n)W(u,n)x(u) ^(t-nTc) 

CDMA decoding of the waveform in FIG. 2 for the transmitter 
is defined in equations (5). Step 1 defines the convolution 
R(n,n-n') of the CDMA pulse waveform with itself in the matched 
filter receiver. Steps 2,3 are the Walsh and PN decoding 
properties- Step 4 uses the matched filter detection theorem to 
derive the estimates of the transmitted symbols. 

CDMA decoding for receiver (5) 
1 Parameters and definitions are defined in 
1,2,3 in equations (4) together with 

R(n,n-n') = convolution of ^(t-nTc) with 

£(t-n'Tc) evaluated at the receiver 
detection times t=nTc 
= R(n) 5(n-n' ) 
= R(n) for n=n' 
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= 0 otherwise 



2 Walsh decoding of channelization codes W(k) 



WW* 



= Nc I 



where 



I = Nc X Nc identify matrix 
W* = conjugate transpose of W 



<W(k) ,W(k' )> = N 5(k-k' ) 



3 PN decoding 
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P(n) P* (n) 



= 1 for all n 



where 



P* = complex conjugate of P 



4 CDMA decoding 



x(k) = P*(n)W*(k,n) z(t)©£(t~nTc) 



It should be obvious to anyone skilled in the 
communications art that these example implementation algorithms 
in equations (1),(2),(3),(4),(5) clearly define the fundamental 
OFDMA and CDMA signal processing relevant to this invention 
20 disclosure and it is obvious that this example is representative 
of the other possible signal processing approaches. 

For cellular applications the encoding algorithms for the 
transmitter describe the implementation of OFDMA and of CDMA 
25 encoding and are the transmission signal processing applicable to 
this invention for both the hub and user terminals, and the 
decoding algorithms for the receiver describes the corresponding 
OFDMA and CDMA receiving signal processing for the hub and user 
terminals for applicability to this invention. 



For optical communications applications the microwave 
processing at the front end of both the transmitter and the 
receiver is replaced by the optical processing which performs the 
complex modulation for the optical laser transmission in the 
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transmitter and which performs the optical laser receiving 
function of the microwave processing to recover the complex 
baseband received signal with the remainder of the signal 
processing functionally the same for the OFDMA and for the CDMA 
5 encoding transmitter and functionally the same as described for 
the OFDMA and CDMA receiving signal processing receiver. 



BRIEF SUMMARY OF THE INVENTION 

10 

This invention introduces the new OWDMA communications 
technology which iir^lements orthogonal multi-resolution complex 
Wavelet division multiple access and is a multi-resolution 
complex Wavelet generalization of OFDMA; introduces the new 
15 application of the multi-scale code division multiple access MS- 
CDMA architecture which integrates MS-CDMA with OFDMA and with 
OWDMA, and introduces the variable gain control over frequency. 

The new OWDMA forms a uniform set of contiguous orthogonal 
20 filters across the available frequency band in one of the 

numerous available architectural options and which implements a 
polyphase filter bank across the available frequency band with 
the basic property that the filters are orthogonal. Each filter 
defines a OWDMA channel for communications. Similar to OFDMA, 
25 the symbol rate within each channel is equal to the channel 

separation 1/NT where N is the number of channels, 1/T is the 
frequency band, and NT is the symbol-to-symbol separation and 
which is made possible by using multi-resolution complex Wavelet 
channelization filters developed in reference [2] . OWDMA filters 
30 are orthogonal in fr.equency which means their frequency spectrums 
are non-overlapping, and they have flat spectrums across each 
channel. 
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OWDMA orthogonality only requires frequency synchronization 
whereas OFDMA orthogonality requires both time and frequency 
synchronization which means that time synchronization errors on 
the user communication channels for the return links to the hub 
5 or access point for cellular communications do not degrade the 
orthogonal separation between the channels as they will for 
OFDMA. 

Sensitivity to frequency synchronization errors is less for 
10 OWDMA primarily because of the non-overlapping of the frequency 
spectrums. The overall tolerance to user-to-user imbalances on 
the return communications channels is better for OWDMA primarily 
because their design keeps the frequency spectrums from 
overlapping in the presence of real operational conditions with 
15 synchronization errors. This allows the forward communications 
link to support a power imbalance between channels to mitigate 
differing ranges and path losses at the user antennas . 

Application of multi-scale code division multiple access 
20 MS-CDMA in reference [2] provides a means to implement the new 
hybrid-Walsh orthogonal CDMA codes in reference [4] over the 
OFDMA/OWDMA channels by spreading the CDMA within each channel 
and over all of the channels such that each user can be spread 
over the complete band. This keeps the chip rate equal to 1/NT 
25 while maintaining the spreading over the fullband 1/T and allows 
the band transmit Tx power to be independently controlled in 
frequency. These two CDMA scales 1/NT, 1/T are generated by MS- 
CDMA in combination with the OFDMA/OWDMA channelization filter 
banks. The 1/T scale is to combat fading and interference similar 
30 to the current CDMA, and the 1/NT scale is for acquisition, 

synchronization, and equalization protection against multi-path 
and provides the flexibility for band power control to provide a 
frequency diversity communications iirprovement . 
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Variable control over the frequency B can be iitiplemented by 
partitioning the CDMA over the channels into separate groups and 
assigning an independent power level to each group of channels. 
OWDMA readily supports differences in power levels between 
5 adjacent channels. Power control is desirable to support 

differences in quality of service SoC, range losses, and path 
loss . 

It should be obvious to one familiar with the CDMA 
10 communications art that the number of scales could be a larger 
number than the two used in this invention disclosure and the 
multi-resolution complex Wavelet design for OWDMA supports the 
partitioning of the frequency band 1/T into several frequency 
scales for the channelization filters, supports individual multi- 
15 resolution complex Wavelet packets for communications in the 

time-frequency domain and which can be integrated into the MS- 
CDMA, and supports separate and segmented communications bands 
simultaneously. 

20 

BRIEF DESCRIPTION OF SEVERAL VIEWS OF THE DRAWINGS 

The above-mentioned and other features, objects, design 
algorithms, and performance advantages of the present invention 
25 will become more apparent from the detailed description set forth 
below when taken in conjunction with the drawings and performance 
data wherein like reference characters and numerals denote like 
elements, and in which: 

30 FIG. 1 is a description of the OFDMA waveform in the 

frequency and time domains over a 20 MHz band for the IEEE 
802. llg standard. 
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FIG- 2 is a description of the CDMA waveform in the 
frequency and time domains over a band B. 

FIG. 3 is a description of the OWDMA waveform in the 
5 frequency and time domains over a band B. 

FIG. 4 is a description of the design requirements on the 
multi-resolution complex Wavelet PSD=|\|/(f) |2. 

10 FIG. 5 is a representative baseband frequency response for 

the multi-resolution complex Wavelet filter and the square-root 
cosine filters for the bandwidth expansion factors =0.22, 0.40. 

FIG. 6 is a representative implementation of MS-CDMA 
15 OFDMA/ OWDMA encoding for transmitter. 

FIG. 7A is a representative MS-CDMA OFDMA/OWDMA transmitter 
implementation block diagram, the. 

20 FIG. 7B is a representative MS -CDMA OFDMA/OWDMA transmitter 

MS -CDMA mapping. 

FIG. 8 is a representative implementation of MS-CDMA 
OFDMA/OWDMA encoding for the receiver. 
25 FIG. 9 is a representative implementation block diagram for 

the MS-CDMA OFDMA/OWDMA receiver. 

DETAILED DESCRIPTION OF THE INVENTION 

30 

This invention introduces the new orthogonal multi- 
resolution complex Wavelet division multiple access OWDMA 
communications technology which is a multi-resolution complex 
Wavelet generalization of OFDMA, introduces the new multi-scale 
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code division multiple access MS-CDMA architecture to integrate 
MS-CDMA with OFDMA, and introduces the new MS-CDMA architecture 
to integrate MS-CDMA with OFDMA, 

S OWDMA uses the noil ti -resolution complex Wavelet waveform 

developed in reference [2] to generate multi-rate orthogonal 
filter banks and waveforms to support communications at the 
critical symbol rate equal to a combined 1/T symbols per second 
for a 1/T Hz frequency band. The critical symbol rate is the 
10 Nyquist sample rate and corresponds to no excess bandwidth a=0 in 
the equation B={H-a)/T which can be rewritten as the bandwidth- 
time product BT=(l+a) that relates filter bandwidth B and the 
symbol rate 1/T for communications supported by this filter. For 
the OFDMA filter in FIG. 1 it is observed BT=l+a=64/52=l . 33 
15 corresponding to a=0-33 since 52 of the 64 DFT filters generated 
are used and for the CDMA filter in FIG. 2 it is observed that 
the bandwidth-time product is about BT=1.25 corresponding to 
a=0.25 for the 3G CDMA cellular communications systems. 

Multi-resolution complex Wavelet design algorithms were 
developed in reference [2] as a means to design polyphase 
multirate filters, quadrature mirror filters (QMF) , perfect 
reconstruction filters, Wavelet iterated filter banks, and 
Wavelet tiling of the time-frequency t-f space. Prior to the 
invention of multi-resolution complex Wavelet design algorithms, 
theorectical studies had not yielded useful realizable filters 
for system applications implementing these architectures as 
summarized by the digital filtering and polyphase research in 
references [5] . [6] and the Wavelet research as summarized in the 
references [7] , [8] , [9] , 10] , [11] . 
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30 



Multi-resolution complex Wavelets defined in equations (6), 
(7), (8) expand the Wavelet analytical formulation to include a 
frequency variable which specifies the center frequency of our 
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new Wavelets and perforins a frequency translation of the real 
mother Wavelet to this center frequency. Currrently, Wavelets are 
real functions of the scale and translation parameters. Multi- 
resolution complex Wavelets are functions of these parameters 
5 plus the frequency variable. This new concept of frequency as an 
additional parameter provides an added degree of flexibility and 
together with the Fourier domain design approach provide an 
entirely new means for deriving these new waveforms as 
generalization of the traditional real Wavelets and their 
10 modification to include the added frequency parameter. With 

frequency translation as the additional parameters the analytical ^ 
formulation of these new waveforms as a function of the baseband 
or mother waveform centered at dc ( dc refers to the origin f =0 of 
the frequency space) is given in equations (6). 

15 

Equations (6) introduce the new multi-resolution complex 
Wavelets. Step 1 in equations (6) is the definition for the 
continuous real Wavelet over the time-frequency t-f space from 
reference [7]^ [9], [12] where v|/(t) is the Wavelet which is a 

20 waveform of finite extent in time t and frequency f over the t-f 
space. Wavelet parameters a,b are the Wavelet dilation and 
translation respectively or equivalently are the scale and shift. 
The \|/(t) without the indices a,b is the mother Wavelet which is a 
real and symmetric localized function in the t-f space used to 

25 generate the doubly indexed Wavelet \|/(t|a,b) where \|/(t|a^b) reads 
the \|/ is a function of time t for the parameters a,b. Scale 

factor |a| has been chosen to keep the norm of the Wavelet 
invariant under the parameter change a,b. Norm is the square root 
of the energy of the Wavelet response. The Wavelets \|/(t|a,b) and 
30 ^f{t) are localized functions in the t-f space which means that 
both their time and frequency lengths are bounded. 
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step 2 in equations (6) defines the digital Wavelet which 
is a Wavelet over the digital t-f space corresponding to 
digitization of the t-f space at the 1/T analog-to-digital A/D 
rate where T is the interval between digital samples and '^i" is 
5 the index over the digital samples. Time "t" for the continuous 
Wavelet in 1 is replaced by the equivalent digital sample number 
"i" corresponding to t=iT at sample i. Wavelets in digital t-f 
space have an orthogonal basis that is obtained by restricting 
the choice of the parameters a,b to the values a=2'P, b=qN2P where 
10 p,q are the new scale and translation parameters and N is 



New multi-resolution complex Wavelets 



(6) 



1 Current continuous Wavelet as a function of the 



mother Wavelet at dc 
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v|/(t|a,b) = |ar'/2 v|/((t-b)/a) 
where a,b are the dilation, shift parameters 



20 



2 Current digital real Wavelet 

Digital Wavelet shift q and scale p 
a - 2P 



b = q N 2^ 



Digital Wavelet as a function of mother Wavelet 
\|/(i|p,q) = 2-P/2 \|/(2-P i-qN) 



where N = number of samples over the Wavelet 
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spacing or repetition interval 
Tg = NT wavelet spacing 
i = digital sampling index 
T = digital sampling interval 
t = iT time index at digital sample i 
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3 Our new digital multi-resolution complex Wavelet 



is a function of the mother Wavelet at dc 



3. a) Generalized complex format 

VdlP/q.k) = \|/(2-Pi-qN) exp [ j 27cf c (p, k) 2'PiT] 



3.b) Multi-rate filter complex format 

vff(i|p,q, k) = 2"P^2 \|/(2-P i-qN)E*(k2P, 2-^1) 

4 Orthogonality equations for our new digital 
complex Wavelet 

Zi \|/(i|p,q,k) \|/*(i|p^q' ,k' ) =N iff p=p',q=q^k=k' 

= 0 otherwise 

5 New digital complex Wavelet for a uniform filter bank 

as a function of the mother Wavelet at dc 
V(i|p=0,q,k) = v|/(i-qN) E*(k,i) 

the spacing or repetition interval Ts=NT of the Wavelets (which 
from a communications viewpoint are symbols) at the same scale p. 
Wavelets at p,q are related to the mother Wavelet by the equation 
in 2 where the mother Wavelet is a real and even function of the 
sample coordinates and which follows directly from 1 for the 
continuous t-f space. 

Steps 3,4 in equations (6) define our new Wavelet in the 
digital t-f space and their orthogonality properties. Our new 
Wavelets are complex generalizations of Wavelets in t-f space 
which enable them to be useful for communications and radar 
applications- This generalization is accomplished 1) by the 
addition of a frequency translation parameter k which controls 
the frequency offset of the Wavelet, 2) by generalizing the 
Wavelet weighted orthonormality properties in step 4 to apply to 
waveforms over the time translation q and over the scales p with 
the inclusion of the frequency translation and where is the 
complex conjugate of \|/ and is required in the orthogonality 
equations since the multi-resolution complex Wavelet \|/ becomes 
complex with the addition of harmonic k, and 3) by their 
characterization and design in the Fourier domain. The frequency 
parameter k controls the frequency translation exp [ j27cfc (p, k) 2" 
"^iT] in 3 for the generalized format and the frequency 
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translation E*(k2P,2"Pi) in 3 for the multi-rate filter format- 
With this frequency translation the analytical formulation 3 of 
these new Wavelets is given as a function of the baseband or 
mother waveform centered at dc corresponding to k=0. Purpose of 
5 the frequency index k is to identify the center frequencies of 
the waveforms at the scale p and time translation q in the t-f 
space. The generalized center frequency fc(P/k) of the frequency 
translated dc waveform at scale p and frequency index reduces to 
fc(p,k)= k2^/NT for application to multi-rate filters. 

10 

Step 5 in equations (6) defines the equations for our new 
Wavelet for a uniform polyphase filter bank which is one of 
several OWDMA candidate architectures- Our new Wavelets in 5 are 
the impulse responses of the corresponding digital symbols 

15 encoded with digital data for transmission which is the synthesis 
filter bank in polyphase theory, and are the filter bank 
detection impulse responses for recovery of the transmitted 
digital symbols by the analysis filter bank in polyphase theory 
assuming matched filter detection. The digital filters are 

20 observed to be the DFT's of the mother wavelets similar to the 
construction of the OWDMA in 3 in equations (2) upon replacing 
the OWDMA pulse waveform £ with the mother Wavelet \|/(i-qN) in 5 
and setting q=0 corresponding to the symbol at t=0. 

25 Equations (7) derive the multi-resolution complex Wavelet 

as a function of the Fourier design coordinates. Design 
algorithms provide a means to design the mother Wavelet in the 
Fourier frequency domain to fit the communications and radar 
specifications. From this mother Wavelet, the Wavelets at the 

30 appropriate scales {p,q, k} are easily found as demonstrated in 
3,5 in equations (6) . Design in the frequency domain means the 
design coordinates specifying the Wavelet are Fourier frequency 
harmonics or coordinates . Step 1 lists the parameters and 
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coordinates and step 2 defines the Fourier harmonic frequency 
design coordinate. 

Step 3 is the DFT representation from equations 3 in 
equations (1) of the real mother Wavelet \|/(i) as a function of 
the DFT Fourier harmonic coefficients {\|/(k')}. Step 4 derives the 
equation for the multi-resolution complex Wavelet as a function 
of the DFT Fourier harmonic coefficients by substituting 3 in 
equations (7) into 3 in equations (6) . 

Digital multi-resolution complex Wavelet is a (7) 
function of the Fourier harmonic design coordinates 

1 Parameters and coordinates 

\|/ = \|f(i) mother Wavelet 

= \|f(i|p=0,q=0,k=0) 
N' = length of the mother Wavelet v|/(i) 
= NL+1 

N = sampling interval of \|/ 

= spacing of \j/ for orthogonality 
L = length of \|/ in units of N 

T = spacing of digital samples 

Tg = N'T length of \|f in seconds 

2 Fourier harmonic frequency design coordinates 

\|i(k' ) = value of the Fourier design harmonic k' /N' T 

= Fourier harmonic design coordinate value 
{k' } = set of hairmonic design coordinates 
for k' = 0, . . . 

corresponding to harmonics {k'/N'T} 

3 Multi-resolution real mother Wavelet definition 

v|/(i) = Zk' v|/(k' ) E*(k' ,i) 
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4 Multi-resolution complex Wavelet definition 
\|/(i|p,q,k) 

= 2"P/2 2k.\|/(k') E*(k'/L+k2P,2-P(i-qN)) 



5 Several fundamental properties follow directly from the 

frequency design approach in 4 in equations (7) . It is 
demonstrated in reference [2] that our multi-resolution complex 
Wavelets are implemented with our design in the Fourier domain 
and our multi-resolution complex Wavelet design remains invariant 
10 under scale changes. It is demonstrated in reference [2] that the 
Fourier frequency domain design in 3 in equations (7) remains 
invariant for all parameter changes and in particular for all 
scale changes. 

15 Multi-resolution complex Wavelet design algorithms are 

illustrated by a representative least-squares LS design for a 
OWDMA polyphase uniform filter bank in FIG. 3- There are two 
categories of LS algorithms and these are the eigenvalue and the 
gradient search that respectively can be reduced to algorithms 

20 equivalent to the original eigenvalue [13] and Remez-exchange 
[14] waveform design algorithms for application to a uniform 
filter bank. 

FIG. 3 is an example of an orthogonal multi-resolution 
25 complex Wavelet division multiple access OWDMA Wavelet polyphase 
filter bank over the band B by setting N in 5 in equations (6) 
equal to the N in 17 in FIG. 3. The band channelization filter is 
a roofing filter h(f) 14 that covers the B frequency band 15 
assigned to OWDMA. Unlike its use in FIG. 1 for the OWDMA with 
30 3=20 MHz this filter is intended to be flat over the frequency 
band B of interest and with a bandwidth 18 equal to 1/T > B. 
Plotted is the power spectral density PSD=|h(f) |^ of this 
channelization filter h(f). The h(f) output is digitized at the 
sample rate 1/T to form the OWDMA multi-resolution complex 
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Wavelet polyphase filter bank 16 which are uniformly spaced 17 
at 1/NT Hz over the 1/T frequency band 1/T. This digital sartple 
rate 1/T is sufficiently large to allow the use OWDMA filters in 
17 in FIG. 3 for communications over the band B 18 with no excess 
5 bandwidth a=0 unlike the OFDMA in FIG. 1 which has a=0.33 and 
CDMA in FIG. 2 with a=0.25. 

This OWDMA polyphase filter bank is ideally decimated which 
means the filter output multi-resolution complex Wavelet sample 
rate 1/Ts is equal to the channel-to-channel spacing 1/Ts=l/NT 
equivalent to stating that there is no excess bandwidth a=0 
within the filter bank. A representative 19 OWDMA multi- 
resolution complex Wavelet \|/ for L=8 is plotted in 20 as a 
function of the time offset expressed in units of the multi- 
resolution complex Wavelet spacing NT. This \|/ was designed by the 
eigenvalue category of LS design algorithms. Our design for this 
topology is immediately applicable to an arbitrary set of multi- 
resolution OWDMA filters through the scaling 3,4 in equations (7) 
which gives the design of our Wavelet at arbitrary scales in 
terms of our design of the mother Wavelet. 

Multi-scale mother Wavelet \\f design for the OWDMA polyphase 
filter bank in FIG. 3 starts by using the frequency design 
template in FIG. 4 to construct the LS error metrics as a 
25 function of the DFT frequency design coordinates {\|/(k' ) in 2 in 
equations (7) which define the mother Wavelet in 3 in equations 
(7). Next minimization search algorithms are developed and used 
to find the values of {\|/(k') that minimize the weighted sum of 
these LS design error metrics equal to the cost function J in 
30 equations (8) for the LS design. Minimizing J with respect to the 
{v|/(k') gives the design values of {\|/{k' } for constructing the 
mother Wavelet in 3 in equations (7) . 



21 



10 



15 



20 



Equations (8) define the LS cost function J constructed 
with the LS metrics with the aid of the frequency design tenplate 
in FIG. 4. Step 1 defines the LS error metrics which are the 
passband error metric stopband error metric ^1(2), 

5 quadrature mirror filter QMF error metric ^(3), intersymbol 
interference error metric |x(4), and the adjacent channel 
interference error metric ^i{5). 

LS cost function J for designing \|/(i) (8) 
10 1 LS metrics 

^(1) = passband LS error metric measures the LS 
error of the passband ripple 24 in FIG. 4 
^1(2) = stopband LS error metric measures the 
stopband atteunation 28 in FIG. 4 
15 ji(3) = quadrature mirror filter QMF LS error 

metric measures the flatness of the sum 
of two contiguous filters over the 
deadband 31 in FIG. 4 
^(4) = orthogonality LS error metric measures the 
20 intersymbol interference ISI between 

overlapping Wavelets of different wymbol 
H(5) = orthogonality LS error metric measures the 
adjacent channel interference ACI from 
nearest neighbor channels 
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2 Metric weighting 

w(n) = metric weight for error metrics n=l,2,3,4,5 in 1 

> 0 

such that Zn w(n) =1 

3 Cost function J 
J = Zn w(n) ^(n) 
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FIG. 4 is the frequency design template for the power 
spectral density PSD=|v|/(f) |^ of the multi-resolution complex 
Wavelet and defines the parameters of interest for the passband 
metric |x(l) and stopband metric in step 1 in equations (8) . 

5 Passband 21 of the wavelet PSD=|\|;(f) |^ is centered at dc (f=0) 
since we are designing the mother Wavelet, and extends over the 
frequency range cOp extending from -(Op/2 to +©p/2 22 in units of 
the radian frequency variable (0=27cfT 23 where T is the digital 
sampling interval defined in FIG. 3. The frequency space extends 

10 over the range of Af=-1/2T to Af=+1/2T which is the frequency 
range in FIG 3 and the mother Wavelet is at the center of the 
frequency band. Quality of the PSD=|\|/(f ) |^ over the passband is 
expressed by the passband ripple 24. Stopband 25 starts at the 
edge 26 of the passbands of the adjacent channels +/-(0a/2 26 and 

15 extends to the edge of the frequency band co=+/-7c 27 

respectively. Stopband attenuation 28 at +/-c)a/2 measures the 
PSD=lY(f) 1^ isolation between the edge of the passband for the 
mother Wavelet and the start of the passband for the adjacent 
Wavelet channel centered at +/-a)s 29. Rolloff 30 of the stopband 

20 is required to mitigate the spillover of the Wavelet channels 

other than the adjacent Wavelet channels, onto the mother Wavelet 
channel. Deadband or transition band 31 is the interval between 
the passbands of contiguous Wavelet channels, and is illustrated 
in FIG. 4 by the interval from (Op/2 to (Da/2 between the mother 

25 Wavelet channel and the adjacent Wavelet channel at ©a. Waveform 
sample rate ©s 32 is the waveform repetition rate. For the LS 
example algorithm, the waveform sample rate is equal to the 
channel-to-channel spacing for zero excess bandwidth. Therefore, 
1/Ts= <i)s/27cT=l/NT which can be solved to give (Os=27c/M for the 

30 radian frequency sampling rate of the filter bank which is 
identical to the Wavelet repetition rate. 
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Equations (8) step 1 QMF LS error metric \i(3) expresses the 
requirement on the deadband that the PSD's from the contiguous 
channels in FIG. 4 add to unity across the deadband 31 [Op , ©3] 
in FIG. 4 in order that the Wavelets be QMF filters. 

5 

Equations (8) step 1 Inter-symbol interference ISI and ACI 
error metrics |x(4),|j,(5) are orthonormality metrics that measure 
how close we are able to design the set of Wavelets to be 
orthonormal over the t-f space, with the closeness given by the 
10 ISI error metric ji(4) and the ACI error metric |i(5). ISI is the 
non-orthogonality error between Wavelets within the same channel 
separated by multiples of the sampling interval 1/MT seconds 
where T is the sample time and M is the interval of contiguous 
samples. Adjacent channel interference ACI is the non- 
15 orthogonality error between between Wavelets within a channel and 
the Wavelets in adjacent Wavelet channels at the same sample time 
and at sample times separated by multiples of the sample 
interval. As observed as noise contributions within each sample 
in a given channel, the ISI is the noise contribution due to the 
20 other received Wavelets at the different timing offsets 

corresponding to multiples of the sanpling interval. Likewise, 
the ACI is the noise contribution due to the other Wavelets in 
adjacent Wavelet channels at the same sampling time and at 
multiples of the sampling interval. 

25 

Equations (8) step 1 ISI and ACI errors are fundamentally 
caused by different mechanisms and therefore have separate 
metrics and weights to specify their relative importance to the 
overall sum of the LS metrics. ISI is a measure of the non- 
30 orthogonality between the stream of Wavelets within a channel as 
per the construction in FIG. 4. On the other hand, ACI is a 
measure of the non-orthogonality between the Wavelets within a 
channel and the other Wavelets in adjacent channels. This means 
the stopband performance metric has a significant iirpact on the 
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ACI due to the sharp rolloff in frequency of the adjacent 
channel, and the ACI metric is then a measure of the residual 
non-orthogonality due to the inability of the stopband rolloff in 
frequency from completely eliminating the ACI errors. 

5 

Equations (8) step 2 defines the weights of the LS error 
metrics when summed to yield the cost function J. These weights 
are real and normalized to sum to unity. They have proven to be 
helpful in the Wavelet design to emphasize the relative 
10 importance of the individual error metric contributions to J. 

Equations (8) step 3 defines the cost function J as the 
weighted sum of the LS error metrics and which is minimized with 
respect to the DFT frequency design harmonics {\|/(k' ) } to select 
15 the best LS choice for the {\|/(k')} to design the mother Wavelet 
in 3 in equations (7) and the channel Wavelets by frequency 
translation in 5 in equations (6) . 

Multi-scale mother Wavelet frequency response in FIG. 5 is 
20 evaluated by irrplementing the LS design algorithms in reference 
[2] that minimize the J in 3 in equations (8) to find the best 
set of DFT frequency design coordinates {^(k')} which give the 
mother Wavelet 3 in equations (7). FIG- 5 plots the PSD 
frequency response for the multi-scale mother Wavelet and the 
25 square-root (sq-rt) raised-cosine (r-c) waveforms with excess 

bandwidth a= 0.22, 0.4 0 which waveformns are extensively used for 
other communications. Plotted are the measured PSD in dB 42 
versus the frequency offset 43 from dc expressed in units of the 
symbol rate, for the new Wavelet waveform 44, the sq-rt r-c 
30 with a=0.22 45, and the sq-rt r-c with a=0.40 46. It is 
believed that the multi-resolution complex Wavelets can be 
designed as a filter with better performance parameters than 
possible with any other known algorithm. 
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OWDMA encoding for the transmitter is defined in 
equations (9) . Step 1 lists parameters and definitions and step 
2 defines the transmitted OWDMA encoded baseband signal z(i) for 
contiguous data blocks and where the. symbol offsets A account for 
5 symbol overlaps over the symbol q data block interval within each 
channel. 



OWDMA encoding for transmitter (9) 

1 Parameters and definitions 

10 h(i) = roofing filter impulse response in time 

for h(f) in 15 in FIG. 3 
\|/(i) = \|/(i |p=0, q=0, k=0) mother Wavelet in 3 
in equations (7) 
= baseband or do k=0 Wavelet at p=0,q=0 
15 N = number of OWDMA filters over the 1/T 

frequency band 
1/T = digital sample rate for OWDMA 

> coitplex Nyquist rate for roofing filter h(f) 
NT = OWDMA Wavelet spacing 
20 1/NT = OWDMA Wavelet output rate 

= OWDMA Wavelet channel separation 
NBT = channels used for data and pilot 

where B is the frequency band in FIG. 3 
N(l-BT) = guard band channels for rolloff 
25 of the h(f) 

d(k,q) = data modulation for user k for data block q 

= encoded amplitude A(k|q) and phase (p(klq) 
x(k|q) = transmitted symbol encoded with d(k) 
= A(k|q) exp(jcp(k|q) ) 
30 Assume the h(f) is flat over the passband for 

both Tx and Rx and can be neglected 

2 Transmitted OWDMA encoded baseband signal z(i) 
Time index field referenced to q=0 
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i = 0, +/-2, . . . 

= digital sampling time index 
io = index over a \|/(i) spacing interval 
= 0,1,2, . . . . ,N-1 
5 A = index over the \|/(i) range in units 

of the \|/(i) spacing interval N 
- (L/2-1) , , - .,-l,0,+l, . . .,+(L/2-l) 
Complex baseband signal 

z(iolq) = complex baseband signal over io 
10 for data block q 

= FWT[ x(klq+A) ] 
= N'^ Sa £k x{k|q+A)\|/(i|q+A, k) 
= N"^ 2a x(k|q+A)v|/(io+AN)E*(k,io) 
z(i) = XqZ(io|q) 
15 where FWT is the fast multi-resolution complex 

Wavelet transform 

3 FWT algorithm for OWDMA encoding in the transmitter 
3. a) FWT pre-calculation FFT"^ 
20 X(io,q+A) = N"^ Sk x(k|q+A) E*{k,io) 

= FFT"^[ N"^ Sk x(k|q+A) E*(k,io) ] 
3.b) FWT post-sum 

z(iolq) = Sa H/(io+AN) >-(iolq+A) 

25 4 Computational complexity of fast algorithm 

Real multiply rate Rm 

R„T = 2 log2(N) + 2 L 
Real add rate Ra 

RaT = 3 log2(N) + 2 L 
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step 3 is the new fast FWT algorithm in this invention 
disclosure for the transmitted OWDMA which consists of the pre- 
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calculation FFT"^ in sub-step 3. a followed by a post-sum in sub- 
step 3.b. 

Step 4 evaluates the real multiply complexity metric RmT and 
5 real add computational complexity metrics RaT in terms of 

multiplies /adds per digital sample for the fast algorithm in step 
3, The first term in these metrics is the complexity of the FFT'^ 
for a base 2 and the second term is the complexity of extending 
the multi-resolution complex Wavelet waveform over L of the 
10 symbol intervals . 

OWDMA decoding for the receiver is defined in equations 
(10) . Step 1 refers to 1,2 in equations (9) for the parameters 
and definitions and defines the OWDMA filtering Wavelet. Step 2 
15 demonstrates Wavelet orthogonality. Estimates of the transmitted 
symbols in step 3 are equal to the FWT"-^ of the received baseband 



signal. 



OWDMA decoding for receiver 

1 Parameters and definitions are defined in 1,2 in 
equations (9) together with 

OWDMA filtering wavelet from 5 in equations (6) 



(10) 
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\|i*(i|q,k) = \|/*(i|p=0,q, k) 
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= \|/(i-qN)E(k, io) 
Assume the h(f) is flat over the passband for 
both Tx and Rx and can be neglected 
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2 Multi-resolution coinplex Wavelet orthogonality from 4 in 
equations (6) 

k , k ' orthogonal i ty 



ZiM/(ilq/K) v|/*(i|q,kM 

= Si \|/2(i-qN)E*(k,i)E(k',i) 

= lio U [lAq v'(io+A)N)] E*(k,io)E(k' ,io) 



= Zio [ 1 ]E*(k,io) E(kMo) 
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= N 8(k-k' ) 
using the completeness property of the 
multi-resolution complex Wavelet from reference [2] 

[2A\|/^(io+ A)N)] E 1 for all io 
5 q,k and q' , k' orthogonality proven in reference 

[2] and restated in 4 in equations (6) 

3 OWDMA decoding derives estimates x (k|q) of x(k|q) for 
data block q from the receiver estimates z (i) of z(i) 

10 x(k|q) = FWT-'[ z(i) ] 

= Si z{i) \|/*(i|q,k) 

4 FWT algorithm for OWDMA decoding in receiver 

4. a) FWT pre-sum 
15 ^(iolq) = 2a z (iolq+A) \|/(io+AN) 

= 2a z (io+AN|q) \|/(io+AN) 
4.b) FWT pre-sum FFT 

x(k|q) = lio X(iolq) E(k,io) 



20 5- Computational complexity of fast algorithm 

Real multiply rate Rm 

RmT = 2 log2(N) + 2 L 
Real add rate Ra 

RaT = 3 log2(N) + 2 L 

25 

Step 4 is the new fast algorithm for the received OWDMA 
which partitions the baseband symbol detection x{k|q) in step 3 

into a pre-sum calculation sub-step 4. a of A.(iolq) followed by a 
sub-step 4.b FFT of this pre-sum, 

30 

Step 5 evaluates the real multiply complexity metric RmT and 
real add computational complexity metrics RaT in terms of 
multiplies/adds per digital sample for the fast algorithm in step 
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3. The first term in these metrics is the complexity of the FFT'^ 
for a base 2 and the second term is the complexity of extending 
the multi-resolution complex Wavelet waveform over L of the 
symbol intervals. 

5 

MS-CDMA parameters and codes are defined in equations (11). 
Step 1 defines the scenario parameters. Step 2 partitions the 
user index u field into the sub-fields uo,Ui of size No,Ni for 
scales 0,1 respectively and which are the indices over the users 
10 within each channel and the indices over the channels within the 
MS-CDMA group and which uniquely represent u as u=Uo+UiNi. 

Step 3 partitions the code chip index n field into the sub- 
fields no^ni of size Nq^Ni for scales 0/1 respectively and which 
15 are the indices over the chips within each channel and the 

indices over the channels of the MS-CDMA group and which uniquely 
represent n as n=no+niNi. 

Step 4 defines the NcXNc MS-CDMA code matrix C whose 
20 elements are C(u,n) where u+1 is the row index and n+1 is the 

column index and where the +1 has been added to correspond to the 
row and column numbering starting with +1. MS-CDMA code vector 

MS-CDMA parameters and codes (11) 
25 1 Scenario parameters 

M = number of communications channels 
of band B 

No = niomber of CDMA chips per channel 
Ni = number of channels in MS-CDMA group 
30 Nc = number of chips in MS-CDMA group = No Ni 

1/To = MS-CDMA chip/symbol rate 

= 1/NT for OFDMA,OWDMA in FIG. 1,2 
x(u, q)= User u in group q 

= x(k|q) when No=l, u=Jc 
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= x(k) when No=l/ u=k, q=0 

User index u for a MS-CDMA group 
u = Index of CDMA users 
= 0,1,..., (Nc"l) 

= (uo,Ui) field representation of u 
= uo + UiNo for fields Uo, Ui 
where the index fields are 
Uo = index of users in a channel in 
the MS -CDMA group 
= 0,1,..., (No-1) 
ui = index of channels in the MS-CDMA group 
= 0,l,...,Ni-l 

Code chip index n for a MS-CDMA group 
n = Index of CDMA code chips 
= 0,1,.,., (Nc~l) 

= (no,ni) field representation of n 
= Ho + niNo for fields no, ni 
where the index fields are 

no = index of code chips within a channel 
in the MS -CDMA group 
= 0,1,..., (No -1) 
ni = index of channels in the MS-CDMA group 
= 0,1,..., ( Ni-1) 

MS-CDMA code matrix C for a MS-CDMA group 
C = Nc X Nc MS-CDMA code matrix 

= [C(u,n)] matrix of elements {C(u,n)} 
= [c(u)] matrix of Ix Nc code vectors 
C(u) = MS-CDMA code vector u which is row u+1 
C counting rows, columns starting with 1 
= [C(u,0), C(u,l), . . . ,C(u, Nc-1)] 
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5 MS-CDMA code 

matrix C construction 

5. a No n- factorable C: algebraic field construction 
C = t C(u,n) ] 
5 = [ C(uo+UiNor no+niNo) ] 

5.b Factorable C: Tensor product (Kronecker) 
construction 
C = Ci®Co 

10 == [Ci(Uuni)Co(uo,no) ] 

where 

® = Kronecker (tensor) product of Ci and Co 
Ci = Ni X Ni MS-CDMA scale 1 code matrix 

= [Ci(ui,ni)] matrix of elements {Ci(Ui,ni)]} 
15 over the user channels within the group 

Co = No X No MS-CDMA scale 0 code matrix 

= [Co(uo/no)] matrix of elements {Ci(ui,ni)]} 
over the user chips within each channel 
in the group 

20 

c(u) for user u is the row u+1 of C- The C is a complex 
orthogonal code matrix - 

25 Step 5 illustrates the construction of MS-CDMA code matrix 

C for non-factorable and factorable C. In 5. a for a non- 
factorable C the algebraic fields for code indices u and chip 
indices n is constructed as the sum of the scale ^^0" algebraic 
fields of the code indices Uq and chip indices no of users within 

30 each channel in the MS-CDMA group, plus the scaled addition of 
the scale "1" algebraic fields of the code indices Ui and chip 
indices ni of users over the channels of the MS-CDMA group 
wherein the scale factor is the size No of the scale "0" algebraic 
fields of indices. For the special case where C is factorable in 
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5.b the C is constructed as is a Kronecker or Tensor product 
C=Ci®Co of Ci and Co where is the Kronecker or tensor product 

and the matrix Ci is the Ni x Ni MS-CDMA scale ^^1" code matrix 
over the user channels within the MS-CDMA group and the matrix Co 
5 is the No x Nq MS-CDMA scale "0" code matrix for the user chips 
within each channel in the MS-CDMA group. 

MS-CDMA representative application to OFDMA in FIG. 1 uses 
the candidate architecture which spreads the CDMA over the M 
10 contiguous 48 data channels or 52 contiguous data plus pilot 
channels and the representative application to OWDMA in FIG. 3 
uses the candidate architecture which spreads the CDMA over the M 
contiguous OFDMA communications channels across B. 

15 MS-CDMA partitions the M channels into channel groups of 

size Ni and provides a scale over the chips within the channels 
and another scale over the channels within this group. Code chip 
length Nc for each user in a channel group is equal to Nc=NoNi 
where Nq is the number of chips in each channel assigned to scale 

20 ^'0" and Ni is the number of channels assigned to scale '^1" 

within the group. Each user has a chip scale "0" and a channel 
scale ^'1". Chip scale ^^0" spreads the data over the chips within 
each channel and channel scale '^1" then spreads the channel chips 
uniformly over the Ni channels with the result that each user 

25 occupies each of the channels within the Ni channel group. 

There could be from 1 to M/Ni channel groups depending on 
the architecture and applications. The use of multiples groups 
M/Ni>l tends to be desirable since the storage requirements and 
30 computational complexity are reduced as the number of groups are 
increased and the spreading advantages within each group tend to 
saturate as the number Ni>16 when the channels within each group 
are spread across the fullband M channels. 
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For No=l there is no CDMA within each channel and the MS- 
CDMA then spreads the signals over each channel within a group 
for both OFDMA and OWDMA to function as a means to spread each 
channel over the fullband M channels and which may be a desirable 
5 architecture when the channels are sufficiently narrow to produce 
a sufficiently long pulse to counter multipath. 

MS-CDMA OFDMA transmitter equations are defined in 
equations (12) for MS-CDMA. Step 1 gives the parameters and 

10 definitions. Step 2 defines the encoding equations for chip no 
for data block q. Sub-step 2. a uses the fast code transform 25 
developed in references [3], [4] to generate the encoded vector. 
Sub-step 2.b uses the inverse FFT'"^ to construct the transmitter 
complex baseband signal z(io|no+qNo) for chip no for data block q 

15 and these signals are combined in sub-step 2.c to generate the 
transmitter signal z(i) for all no,q. 

Step 3 evaluates the real multiply complexity metric RmT 
and real add computational complexity metrics RaT in terms of 
20 multiplies/adds per digital sample for the fast algorithm in step 
3. The first term in these metrics is the complexity of the FFT"^ 



MS-CDMA OFDMA transmitter equations 



(12) 



1 Parameters and definitions are defined in 1 in 
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equations ( 2), (9), (11) together with 
MS -CDMA group index field {g} 



g = 0,1,. . . , (M/Ni-1) where M/Ni = number 



of MS-CDMA groups over the M data channels 
Frequency harmonic k 
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k = k(ni|g,q) 
MS-CDMA chip index fields 



n = no + niNo 



= 0,1,2 



,Nc-l 
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no = 0, 1, 2 



= 0,1,2 



,No-l scale 0 chips in channel 
,Ni-l scale 1 chips over channels 
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Nc = NoNi number of chips in a MS-CDMA group 
MS -CDMA channel code index ni 

ni - ni(k|g,q) function of k for a given g,q 
User symbol x 

X = x(u|g,q) 

= A(ulg^q) exp j(p(u|g,q) 
MS-CDMS code orthogonality 

CC* = Nc I 

PN spreading code P(n) for chip n in group g 
P(n|g,q) = exp ( j(p(n|g,q) 

Assume the band and pulse filtering can be 
neglected on both Tx and Rx which is 
equivalent to the assumption £(i)=l for all i 

MS -CDMA OFDMA encoding for complex baseband signal z(i) 
2. a) FCT MS-CDMA encoding for data block q 
Z(n|g,q) FCT [ x(u|g,q) ] 

= Nc""^ Su x(u|g,q)C(u,n)P(n|g,q) 
where FCT is the fast code transform 
2.b) FFT"^ OFDMA encoding for chip no of data block q 
z(io|no+qNo) = FFT"" [ Z(n|g,q) ] 

= Zk Z(n|g,q)E*(k,io) 
where n = no + niNo 
ni = ni (k| g,q) 
Z(n|g,q)=Z (no+ni ( k 1 g, q) No) 
2.c) Complex baseband transmitter signal z(i) 
z(i) = 2(o) z(io|no+qNo) for (o) = (no+qNo) N 



3 Computational complexity 
Real multiply rate Rm 

RmT = 2 log2(N) 
Real add rate Ra 

RaT = 3 log2(N) + 2 log2(Nc) + 2 
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for a base 2 and the second term is the complexity of the FCT 
assuming the FCT does not require any multiplies. 

MS-CDMA OWDMA transmitter equations are defined in 
5 equations (13) for MS-CDMA. Step 1 lists the parameters and 
definitions. Step 2 defines the encoding equations for chip no 
for data block q. The FCT on the symbols in sub-step 2. a yields 
the encoded data block Z(n|g,q) and the FWT on this output in 
sub-step 2.b yields the transmitter complex baseband signal 
10 z(iolno+qNo) for chip no for data block q and these signals are 

combined in sub-step 2.c to generate the transmitter signal z(i) 
for all no,q. 

MS-CDMA OWDMA transmitter equations (13) 
15 1 Parameters and definitions are defined in 1 in 

equations (9), (12) together with the range of q 
q = 0 for the data block being addressed 
= +/-1 contiguous data blocks 
Assume the h(f) is flat over the passband for 
20 both Tx and Rx and can be neglected 

2 MS-CDMA OWDMA encoding for chip no for data block q 
2. a) FCT MS-CDMA encoding for data block q 
Z(n|g,q) = FCT[ x(u|g,q) ] 
25 = Nc"' Xu x(u|g,q)C(u,n)P(n|g) 

2.b) FWT OWDMA encoding for chip no of data block q 
z(iolno+qNo) = Esq 2a FWT [ Z (n+A| g, q+5q) ] 
= N"^Ek2:5qZA Z(n+A|g,q+5q)v|/{ilno+A,k) 
= N'^ZkEsqlA Z(n+A|g,q+6q)\)/(io+AN)E*(k,io) 
30 where the overlap conditions on 5q are 

5q = A for No=l 

= sign (A+no)L I A+no I /NoJ otherwise 
sign(o) = numerical sign of (o) 
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L(o)J = rounded down integer value of (o) 
where the overlap index A = 0, +/-1, . 
. . .,+/-(L/2-l) 

2. c) Complex baseband transmitter signal z(i) 

5 z(i) = E(o) 2(iolno+qNo) for (o) = (no+qNo)N 

3 FWT algorithm for OWDMA encoding in the transmitter 

3. a) FWT pre-calculation FFT"^ 

X(io|no+A,q) 

10 = N"^ ZkZ(n+A|g,q)E*(k,io)] 

= FFT^'[N-' ZkZ(n+A|g,q)E*(k,io)]] 
3.b) FWT post-sxjm 

z(io|no) = £q La Y(i- (no+A)N)X(iolno+A,q) 

15 4 Computational complexity for fast algorithms in 2. a, 3 

Real multiply rate Rm 

RmT = 2 log2(N) + 2 L 
Real add rate Ra 

RaT = 3 log2(N) + 2 L + 2 logalNc) + 2 

20 

step 3 is the new fast FWT algorithm in this invention 
disclosure for the transmitted OWDMA which consists of the pre- 
calculation FFT"^ in siab-step 3. a followed a post-sum in sub-step 
3.b of the product from sub-step 3,a with the corresponding 
25 Wavelet overlaps over the q data block interval. 

Step 4 evaluates the real multiply complexity metric RmT and 
real add computational complexity metrics RaT in terms of 
multiplies/adds per digital sample for the fast algorithms FCT 
30 and FWT in steps 2,3. The first term in these metrics is the 

complexity of the FFT"^ for a base 2, the second term 2L is the 
complexity of extending the multi-resolution complex Wavelet 
waveform over L of the MT symbol intervals, and the remaining 



37 



terms are the complexity of the FCT assuming the FCT does not 
require any multiplies. 

MS-CDMA OFDMA/OWDMA encoding for the transmitter in FIG. 6 
5 is a representative implementation of the MS~CDMA OFDMA encoding 
algorithms in equations (12) and the MS-CDMA OWDMA encoding 
algorithms in equations (13) . Signal processing starts with the 
input stream of data encoded symbols x(u|g,q) 40 from the 
transmitter symbol encoder 52 in FIG. 7A and defined in 1 in 
10 equations (12) for both OFDMA and OWDMA - 

FIG. 6 MS-CDMA encoding 41 implements the fast code 
tramsform FCT encoding defined in sub-step 2. a in equations 
(12), (13) for the MS-CDMA encoding and PN cover or spreading 
15 encoding to generate the Nc outputs x (u | g, q) C (u, n) P (n | g) for each 
MS-CDMA group g and data block q. For each group these outputs 
are siommed 42 over u to generate the encoded vector Z(n|g,q) in 
43. 

20 OFDMA processing 44 performs an FFT'^ on the received set of 

vectors Z(n|g,q) and a summation to implement sub-steps 2.b,2.c 
in equations (12) and the output is band filtered 46 to generate 
the MS-CDMA OWDMA encoded complex baseband signal z(i) in 47. 

25 OWDMA processing 45 performs an FWT on the received set of 

vectors Z(n|g,q) and a summation to implement sub-steps 2.b,2.c 
in equations (13) and the output is band filtered 46 to generate 
the MS-CDMA OWDMA encoded complex baseband signal z(i) in 47. 

30 Outputs z(io|no) 47 from the MS-CDMA OFDMA and MS-CDMA OWDMA 

are digital-to-analog DAC converted 48 and handed off to the 
analog front end 49 as the complex baseband analog signal z(t) in 
49. 
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MS-CDMA OFDMA/OWDMA transmitter description in FIG. 7 
presents a block diagram in FIG. 7A and a representative MS-CDMA 
mapping in FIG. 7B. FIG. 7A is a representative transmitter 
implementation of the MS-CDMA OFDMA and MS-CDMA OWDMA encoding in 
5 FIG. 6. The transmitter block diagram in FIG. 7A includes the 
FIG. 6 MS-CDMA OFDMA/OWDMA encoding in an abbreviated format 54 
in FIG. 7A. FIG. 7A signal processing starts with the stream of 
user input data words. Frame processor 51 accepts these data 
words and performs the encoding and frame formatting wherein CRC 
10 is a cyclic redundant code for error detection, and passes the 
outputs to the symbol encoder 52 which encodes the frame symbols 
into amplitude (Ampl.) and phase coded symbols x(u|g,q) 53 which 
are the input to the MS-CDMA encoding 55 and which is 41,42 in 
FIG. 6. 

15 

MS-CDMA FCT encoding outputs Z(n|g^q) 56 are handed over to 
the OFDMA and OWDMA processing 57 which performs an inverse FFT^^ 
followed by a band filtering for OFDMA which is 44,46 in FIG. 4 
and performs an inverse FWT'^ followed by a band filtering which 
20 is 45,46 in FIG. 6. This complex baseband signal z(i) in 47 in 
FIG. 6 is digital-to-analog DAC converted 59 and the output 
complex baseband analog signal z(t) 60 is handed off to the 
analog front end 61. 

25 The z(t) is single sideband upconverted, amplified, and 

transmitted (Tx) by the analog front end 61 as the real 
waveform v(t) 62 at the carrier frequency fo whose amplitude is 
the real part of the complex envelope of the baseband waveform 
z(t) and the phase angle <|> accounts for the phase change from the 

30 baseband signal to the transmitted signal. Output waveform 62 

from the analog front end is the Tx waveform from the Tx antenna. 

FIG. 7B illustrates a representative MS-CDMA uniform 
mapping of each data symbol over frequency, time, antennas, and 
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beams of a cellular conimunications transmitter. Multiple 
antennas and beams are used when a multiple-input -multiple-output 
MIMO communications link is being implemented. The algebraic 
field construction of the algebraic index fields for the codes 
5 and chips for a 2-scale MS-CDMA construction of a non-factorable 
code matrix C in equations (11) represented by 151 and 152 in 
FIG. 7B for the algebraic chip indices no, ni is continued in FIG. 
7B to include the algebraic chip index fields for chips nz in 
153 over the frequency bands, chips in 154 over the data blocks, 
10 chips n4 in 155 over the transmit antenna beams, and chips ns in 
156 over the transmit antennas. The corresponding algebraic code 
indices are respectively Uo, Ui, U2, U3, U4, U5 and the MS -CDMA 
code length is Nc=NoNiN2N3N4N5 chips . 

It should be obvious to anyone skilled in the 
communications art that this example implementation in FIG, 6,7 
clearly defines the fundamental MS-CDMA OFDMA and MS-CDMA OWDMA 
signal processing relevant to this invention disclosure and it is 
obvious that this example is representative of the other possible 
signal processing approaches. 

MS-CDMA OFDMA receiver equations are defined in 
equations (14) . Step 1 lists the parameters and definitions and 
the assumption that the band and pulse filtering can be 
neglected. Step 2 defines decoding of the received chip signal to 
derive the estimate x(u|g,q) of the transmitted symbol x(u|g,q). 
Sub-step 2. a derives the estimate for the encoded symbols using 
the FFT and sub-step 2.b uses the FCT"^ on this estimate to derive 
the transmitted symbol estimate. 

MS-CDMA OFDMA receiver equations (14) 
1 Parameters and definitions are defined in 1 in 
equations (2), (9), (11) together with 

2 (i) = receiver band filter estimate of 
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the transmitted complex baseband signal 2(i) 
Assume the band and pulse filtering can be 
neglected on both Tx and Rx which is 
equivalent to the assumption 2(i)=l for all i 

5 

2 MS-CDMA OFDMA decoding to derive data symbol x (u|g,q) 
2. a) FFT of each received chip vector 2(io|no+qNo) 
Z (n|g,q) = estimate of transmitted Z(n|g,q) 
= FFT[ z (io|no+qNo) ] 
10 = lio z (iolno+qNo) E(k,io) 

where z(io|no+qNo) = z (i=io+ (no+qNo) N ) 
2.b) FCT"^ of each group g encoded vector Z(n|g,q) 

x(u|g,q) = FCT-'[ Z (n|g,q=0) ] 

= Sn Z (n|g,q)P*(n|g)C*(u,n) 

15 

3 Computational complexity is the same as calculated 
in 3 in equations (12) for the transmitter 

MS-CDMA OWDMA receiver equations are defined in equations 
20 (15) . Step 1 lists the parameters and definitions and the 

assumption that the band and pulse filtering can be neglected, 

MS-CDMA OWDMA receiver equations (15) 
1 Parameters and definitions are defined in 1 in 
25 equations (2) , (9) , (11) , (14) together with 

z (i) = receiver band filter estimate of 

the transmitted complex baseband signal z(i) 
Assume the h(f) is flat over the passband for 
both Tx and Rx and can be neglected 
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2 MS-CDMA OWDMA decoding to derive data symbol x (u|g,q) 
2. a) FWT'^ of each received signal z (i) 

Z(n|g,q) = estimate of transmitted Z(n|g.q) 
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== FWT'^[ z (i) ] 

= Zi z (i)\|/*(ilq,k) 

= EioEa z (io|no+(q+A)N)v(io+AN)E(k,io) 
where z(iolno+qNo) = z (i=io+ (no+qNo)N ) 
k g which means k specifies g 
no,k n=no+ni {k|g,q) 



2.b) FCT ^ of each group g encoded vector Z(n|g,q) 
k (u|g,q) « FCT-U Z (n|g,q) ] 
10 « Sn Z (n|g,q)P*(n|g)C*(u,n) 



3 FWT'^ algorithm 
2. a) FWT"^ pre-sum 

A.(iolg,q) = 2a z (io|no+(q+A)N)\|f(i-(no+A)N) 
15 where Hq+A g using the boundary conditions 

in 2.b in equations (13) 
2.b) FWT"^ pre-sum FFT 

Z(n|g,q) = Zio FFT [ A,(iolg,q) ] 
= 2io X(iolgrq)E(k,io) 

20 

3 Computational complexity is the same as calculated 
in 3 in equations (13) for the transmitter 

Step 2 defines decoding of the received chip signal to derive 
25 the estimate x (u|g,q) of the transmitted symbol x(u|g,q). Sub- 
step 2. a derives the estimate for the encoded symbols using the 
FWT"^ and sub-step 2.b uses the FCT"^ on this estimate to derive 
the transmitted symbol estimate. Step 3 is the new fast FWT 
algorithm in this invention disclosure for the received MS-CDMA 
30 OWDMA which consists of the pre-sum in sub-step 3. a followed an 
FFT on this pre-sum. 
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MS-CDMA OFDMA/OWDMA decoding for the receiver in FIG. 8 is 
a representative iir^jlementation of the MS-CDMA OFDMA decoding 
algorithms in equations (14) and the MS -CDMA OWDMA decoding 
algorithms in equations (15) . Signal processing starts with the 
5 input intermediate frequency IF signal after being single 

sideband downconverted and synchronized in frequency 63 from the 
receiver front end 72 in FIG. 9. This input signal is band 
filtered 64 and handed off to the analog-to-digital converter ADC 
or A/D 65 whose digital output is the received estimate z(i) 66 

10 of the transmitted complex baseband signal z(i) 60 in FIG. 7. 
For OFDMA 67 this received signal is processed by an FFT to 
derive estimates Z(n|g,q) 68 of the MS-CDMA encoded signal 
Z(n|g,q) 56 in FIG. 7. For OWDMA 67 this received signal is 
processed by an inverse multi-resolution coitplex Wavelet 

15 transform FWT"^ to derive estimates Z(n|g,q) 68 of the MS-CDMA 
encoded signal Z(n|g,q) 56 in FIG. 7 implementing the fast 
algorithm defined in 2 in equations (15) . Recovered estimates 
Z(n|g,q) 68 are processed by the inverse fast code transform 
FCT"^ 69 to derive estimates x(u|g,q) 70 of the transmitted data 

20 symbols x(u|g,q) 53 in FIG. 7 for hand over to the symbol 
decoder. 

MS-CDMA OFDMA/OWDMA receiver block diagram in FIG. 9 is a 
representative receiver implementation of the MS-CDMA OFDMA and 

25 MS-CDMA OWDMA decoding in FIG. 8. The receiver block diagram in 
FIG. 9 includes the FIG. 8 MS-CDMA OFDMA/OWDMA decoding in an 
abbreviated format 74 in FIG. 9. FIG. 9 signal processing starts 
with the received Elx waveform 71 from the transmitter 62 in FIG. 
7. Received (Rx) signal v(t) 71 is an estimate of the transmitted 

30 signal v(t) 62 in FIG. 7 received with errors in time At, 
frequency Af, phase A9, and with an estimate z (i) of the 
transmitted complex baseband signal z(t) 60 in FIG. 7. This 
received signal v(t) is amplified and downconverted by the analog 
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front end 72 and single side band SSB downconverted and 
synchronized 73 and handed over to the digital-to-analog 
conversion DAC processing 75 for band filtering by h(f) and 
digitization to generate the baseband signal z(i) 76 which is the 
5 received estimate of the transmitted signal z(i) 60 in FIG. 7. 
Timing synchronization could be implemented in the DAC. 

Outputs z (i) are processed by the MS-CDMA OFDMA/OWDMA 
decoding to derive estimates x(u|g,q) 79 of the transmitted 
symbols x(u}g,q) 53 in FIG. 7 and part of the information is 
handed off to the synchronization (sync) processor. For the 
inverse MS-CDMA'^ OFDMA*"^ the processing 79 consists of the FFT to 

recover estimates Z (n|g,q) of Z(n|g,q) followed by an inverse 
FCT"^ to recover x(u|g,q). For the inverse MS-CDMA"^ OWDMA'^ the 
processing 79 consists of the inverse FWT"^ to recover estimates 

Z(n|g,q) of Z(n|g,q) followed by an inverse FCT"^ to recover 
X (u|g,q) and the fast transform in 2 in equations (15). Outputs 
X (u|g,q) are processed by the symbol decoder 80 and the frame 
processor 81 for handoff as the received Rx data 82. 

It should be obvious to anyone skilled in the 
communications art that this example implementation in FIG. 8^9 
clearly defines the fundamental MS-CDMA OFDMA and MS-CDMA OWDMA 
signal processing relevant to this invention disclosure and it 
is obvious that this example is representative of the other 
possible signal processing approaches. 

Variable power control across the frequency band can be 
implemented by assigning each group g of transmit Tx signals 
30 their own power level P(g). Each MS-CDMA group g occupies a 

subset of the channels over the frequency band B consisting of 
Ni channels which means that the users within group g are 
transmitted with the same Tx power. On receive each group g of 
channels is processed separately so there is no cross-talk 
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between the users in the different groups. The OWDMA was designed 
to support large dynamic range imbalances between channels which 
could be present with power level control. MS-CDMA OFDMA/ OWDMA 
variable power control over the frequency subbands corresponding 
5 to the MS-CDMA groups supports the potential for diversity 

improvements by allocation of the available power to emphasize 
the *best' set of available subbands, which subbands are not 
required to be contiguous as well as supporting the simultaneous 
support of multiple users with differing power requirements due 
10 to range, multi-path, and path attenuation effects. 

A second configuration for variable power control is 
described in reference [2] and increases the flexibility of power 
control to all of the individual channels. 

15 

Preferred embodiments in the previous description is 
provided to enable any person skilled in the art to make or use 
the present invention. The various modifications to these 
embodiments will be readily apparent to those skilled in the art, 
20 and the generic principles defined herein may be applied to other 
embodiments without the use of the inventive faculty. Thus, the 
present invention is not intended to be limited to the 
embodiments shown herein but is not to be accorded the wider 
scope consistent with the principles and novel features disclosed 

25 herein. 



30 



45 



